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ABSTRACT 


A  simple,  empirical  method  has  been  developed  for  pre¬ 
dicting  at  supersonic  speeds  the  normal  force  coefficient,  C^, 

(including  carryover)  of  very-low-aspect  ratio  lifting  surfaces 
mounted  on  bodies  of  revolution.  Predicted  values  of  C^  using 

this  method  are  shown  to  be  in  good  agreement  with  test  data 
obtained  on  both  thick  and  thin  surfaces,  at  Mach  numbers  from 
about  2.5  to  7.7  and  angles  of  attack  to  24°. 

(.  -  <" , 

SYMBOLS  AND  NOMENCLATURE 


A  ,A  ,A 
c’  I’  o 


A.R. 


cross-sectional  areas  of  the  forebodies  of  the 
inlets,  the  inlets,  and  the  freestream  tube 
captured  by  the  inlets  of  ramjet  missiles,  re¬ 
spectively 

2  2 

aspect  ratio  =  b  /S^,  or  b  /Sw 


' ! 


in 


2 


b/2 


exposed  semi-span  of  a  lifting  surface  mounted 
on  a  body  of  revolution 


cross-flow  drag  coefficient 


normal  force  coefficient,  normal  force/qS 

acN/aa 


-  CM  at  0  =  0° 


root  chord 


in 


per  rad 


in 
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d 

E 


“B»*W 

M 

q 

s 

SE’SW 

t 

X 

X 

c.p. 

a 

0 

7 

0 

Subscripts 

B 

BE 

BW 

E 

I 

W 


reference  diameter;  diameter  of  body  on 
which  lifting  surfaces  are  mounted 

complete  elliptic  integral  of  second  kind  with 

2  2  4 

modulus  (1-0  cot  -A.  )* 

Morikawa's  interference  factors 

Mach  number 
dynamic  pressure 
reference  area,  nd  / 4 

total  planform  area  of  housings  (wings) 

average  thickness  of  lifting  surface 

body  station;  X  =  0  at  nose  tip  of  body 

center-of-pressure  location 

angle  of  attack;  angle  between  the  velocity 
vector  and  the  longitudinal  axis  of  the  body 


ratio  of  specific  heats;  y  =  1.4  used  herein 

leading  edge  sweep  angle  for  delta  wings 

aerodynamic  roll  angle;  at  0  =  0°  tne  lifting 
surfaces  are  normal  to  the  plane  of  a 


in 


lbs /in 
in2 


in 

in 

in 

in 

deg 


2 


deg 

deg 


body  alone 

body-housing  combination 
body-wing  combination 
housing 

refers  to  inlet  forebody  and  internal  lift  as  in  AC^ 
wing 
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INTRODUCTION 


The  requirement  for  compactness  in  U.  S.  Navy  missile  designs  results 
frequently  in  configurations  which  incorporate  thick  lifting  surfaces  of 
very-low-aspect  ratio.  These  surfaces  are  invariably  thick,  e.g.,  Figures  1 
and  2,  because  they  are  used  to  house  electronics  and  hydraulics  or  serve  as 
ducts,  as  in  the  case  of  side-mounted  inlets  on  ramjet  missiles.  Current 
requirements  on  missile  speed  have  increased  to  regions  where  guidance  for 
making  aerodynamic  estimates  for  these  surfaces  is  not  available,  either  from 
theory  or  experiment. 


Empirical  estimation  of  the  normal  force  coefficient,  C^,  and  center-of- 


pressure  location, 


X 

c.p. 


) 


for  these  surfaces  is  difficult  because  the  shapes 


are  usually  unique  for  each  new  missile  design  and,  therefore,  the  limited 
test  data  available  are  invariably  for  shapes  that  are  quite  different  from 
the  proposed  shape  in  a  new  missile  design.  Existing  empirical  methods1 
have  been  derived  for  a  specific  class  of  surfaces  and  apply  to  the  lower  end 
of  the  Mach  number  range  of  interest  in  this  presentation. 


Simple  theoretical  methods  that  have  been  used  (with  limited  success)  do 
not  take  into  account  the  effects  of  Mach  number.  Those  that  do,  are  not 
applicable  at  the  very  low  values  of  aspect  ratio  inherent  to  these  types  of 
surfaces.  These  concerns  have  been  expressed  for  some  time.2’3 


A  need  exists,  therefore,  for  either  an  empirical  data  base  for  a  more 
general  class  of  low-aspect-ratio  lifting  surfaces  or  a  simple  predictive 
method  that  is  adequate  in  preliminary  design  for  predicting  and  of 

this  class  of  surfaces  in  speed  ranges  from  moderate  supersonic  to  hypersonic. 

A  simple,  empirical  predictive  method  for  estimating  CN  for  very-l'ow- 

aspect  ratio  surfaces  is  presented  herein.  It  is  shown  to  provide  estimates 
of  CN  that  are  adequate  for  preliminary  design  for  a  variety  of  thicknesses 

and  shapes  and  a  wide  range  of  Mach  numbers  (M  =b  2.5  to  7.7)  and  angles  of 
attack  (a  to  24°). 


OBJECTIVE 

The  objective  of  this  study  was  to  determine  a  simple  method  for  esti¬ 
mating  in  preliminary  design  the  aerodynamic  normal  force  ^icient  of 

very-low-aspect  ratio  lifting  surfaces  (and  body-wing  car  ,  at  moderate 

supersonic  to  hypersonic  speeds  and  to  moderate  angles -of-atti.  . 


METHODOLOGY 


A.  BACKGROUND 

The  method  presented  herein  for  predicting  of  very-low-aspect-ratio 

lifting  surfaces  evolved  from  observations  of  the  experimental  lifting  char¬ 
acteristics  of  thick  surfaces  such  as  those  depicted  by  the  housings  on  the 
wind  tunnel  model  shown  in  Figure  1.  This  model  is  representative  of  an 
Integral  Rocket-Ramjet  (IRR)  missile.  It  was  tested  by  APL/JHU  in  order  to 
compile  aerodynamic  design  information  for  components  of  this  class  of  con¬ 
figurations  since  empirical  methods  for  predicting  and  for  this 

type  of  configuration  and  combinations  of  components  were  not  available. 

Hart's  empirical  curves1  had  been  shown  to  provide  good  predictions  for  low- 
aspect-ratio  surfaces  at  3.0,  but  these  curves  had  been  derived  mostly 

for  wings  that  were  primarily  thin  surfaces.  The  applicability  of  this 
method  to  thick  surfaces  and  to  higher  Mach  numbers  was  therefore  not  known. 

Simple  theoretical  methods  that  account  for  the  Mach  number  variation  of 
AC^  noted  from  test  data  are  not  applicable  at  the  very-low-aspect  ratios  of 

interest  herein;  those  derived  for  aspect  ratios  approaching  zero  do  not 
account  for  the  Mach  number  effects.  This  is  demonstrated  in  Figure  3  wherein 
the  test  values  of  AC  of  the  IRR  E  housings  (AC  minus  internal  momentum) 

\  1  E 

are  compared  with  two  simple  theoretical  methods,  viz:  modified  Newtonian 
theory,  plus  wing-body  carryover,  i.e., 

"Vv4t  (l  -yh  +  rsi"2°  <» 

and  slender  wing  theory6  plus  cross  flow  as  recommended  by  Flax  and  Lawrence,3 
i.e. , 


TT  A.R. 

2  X  57.3 


a  +  Cn  sin  a 


(KW  +  V  s' 


(2) 


The  value  of  C^  =  1.0  was  used  in  these  calculations  following  the  recom- 
c 

mendation  of  Flax  and  Lawrence  for  the  case  of  rounded  tips.  This  number, 
however,  could  be  something  other  than  1.0  according  to  Hoerner.6  The 
Morikawa  carryover  factors'  were  used  in  Equations  1  and  2,  and  are  used 
throughout  this  analysis. 

B.  APPROACH 

Test  data  obtained  on  both  thick  and  thin  wings  in  various  APL/JHU  aero¬ 
dynamic  research  and  exploratory  development  programs  were  the  primary  source 
of  data  for  the  development  of  the  empirical  method  presented  herein.  Selected 
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NASA  data  were  also  used.  Sketches  of  the  housing  and  wing  configurations 
used  in  the  analysis  are  given  in  Figures  4,  5,  and  6;  the  sources  for  the 
test  data8*16  are  noted  in  the  figures  for  each  configuration. 


In  all  cases,  che  wing  (or  housing)  data  were  obtained  from  tests  con¬ 
ducted  with  cruciform  wing-body  and  with  body  alone  configurations.  The 
wing-body  was  roll  oriented  at  0  =  0°,  i.e.,  one  pair  of  wings  in  the  angle- 
of -attack  plane.  The  test  data  then  are  derived  from  AC.,  =  CN  -  CN  and 

W  BW  NB 

thus  wing-body  carryover  is  included  in  the  wing  lift. 

The  general  approach  in  deriving  and  evaluating  the  present  method  using 
the  test  data  discussed  above,  is: 

1.  Values  of  P  C„  were  extracted  from  test  data  obtained  on 

“w 

the  wing  configurations  sketched  in  Figures  4,  5,  and  6. 


2.  Correlation  curves  of  p  C.T  were  deduced  from  the  test  values 

N 

“w 

as  follows:  p  =  F  (p  A.R.)  for  rectangular  wings, 


P  Cjj  =  F  (p  cot«/V  )  for  delta  wings,  and 


p  CN  *  constant  for  thick  wings. 

“w 

A  comparison  of  the  derived  curves  with  appropriate  linear  and  slender  wing 
theories  is  given. 

3.  These  empirically  derived  curves  were  then  used  to  calculate  the 
values  of  AC„  for  the  29  Mach  number-configurational  conbinations  used  in 

the  analysis.  Comparisons  with  test  data  are  given  to  demonstrate  the  ade¬ 
quacy  of  the  present  method. 


RESULTS 

A.  PROCEDURE  FOR  EXTRACTING  P  C„  FROM  TEST  DATA 

N 

“W 

Values  of  P  C^  that  provide  a  good  representation  of  the  test  data  in 
°W 

the  range  of  angle  of  attack  tested  were  derived  by  first  linearizing  AC.,  vs  a 
as  demonstrated  in  Figure  7  and  then  extracting  C^  from  the  linearized  values 

aw 
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-H 


I 


of  AC„  as 

”w 


follows : 


(57.3)  (S/Sw) 
*W  +  KB 


per  rad. 


(3) 


where  ACM  =  CN  -  C  and  this  includes  mutual  body -wing  carryover.  The 
\  BW  NB 

carryover  factors  Ky  and  were  obtained  from  Morikawa's  charts.  Reference  7; 
Morikawa's  values  of  Ky  for  rectangular  wings  were  used  for  the  configura¬ 
tions  that  are  nearly  rectangular.  In  the  linearization  of  AC  vs  a,  more 

emphasis  was  given  to  obtaining  a  representation  of  AC  at  the  moderate  to 

“w 

higher  values  of  a  than  at  the  lower  values  according  to  the  objective  of 
this  investigation. 


The  E^,  E^,  and  E^  configurations  of  Figure  4  have  flow  through  the 

inlet-duct  system  and  thus  AC^  for  these  configurations  include  internal 

£ 

lift.  The  lift  attributed  to  the  inlet  forebody  and  internal  momentum  was 
subtracted  from  the  total  lift  of  these  housings  in  order  to  obtain  AC 

since  we  are  only  interested  in  the  external  lift.  Thus,  for  these  con¬ 
figurations, 


A  value  of  A  /A  =  1.0  was  used  in  these  calculations  since  the  internal  shock 

was  not  expelled  for  the  cases  considered.  A,  and  A  are  the  combined  cross 

I  c 

sectional  areas  of  the  inlets  and  inlet  forebodies,  respectively. 


Finally,  the  derived  slopes  were  expressed  in  the  usual  functional  forms 
found  in  design  charts,  i.e., 

P  CN  =  F  (p  A. R. ) 

°W 

for  rectangular  wings,  and 


P  C^  *  F  (p  cot  -A.  ) 

“w 


for  delta  wings. 
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B.  CORRELATION  CURVES  OF  B  C„ 

N 


1.  Nearly-Rectangular  Wings 

The  "best  fit"  values  of  {3  C^  deduced  from  the  test  data  on  the  nearly- 

°W 

rectangular  housings  and  wings  of  Figures  4  and  5  are  plotted  in  Figure  8  as 
a  function  of  1/p  A.R.  For  comparison,  the  values  of  (3  C^  predicted  from 


a, 


W 

linear  and  slender  wing  theories.  Reference  5,  for  rectangular  wings,  are 
also  shown  in  Figure  8,  i.e., 


P  CN  ~  4  (1  -  2  p  k  R 
aW 


p  A.R.  >  1 


P  CN  ~  n  [<2  P  A.R.  )  8in  1 


“w 


+  (1  + 


i— 


P  A.R. 


P  A.R.  +  (p  A.R.  -  2)  cosh 


1  1 


)  V  1  -  (P  A.R. ) 


] 


P  A.R. 


i  <  P  A.R.  <  1 


and. 


P  CN  =  -  (p  A.R.) 

“w 


P  A.R.  <  1/2 
(Slender  Wing) 


It  is  seen,  from  Figure  8,  that  the  difference  between  experiment  and  theory 
(given  by  these  simple  methods)  is  very  large  for  P  A.R.  <■  0.67 
[(1/P  A.R.)  «  1.5]. 

Note  that  the  theoretical  values  of  C^  are  lift  curve  slopes  at  o  =  0° 

a 

whereas  the  test  values  are  the  mean  values  of  C„/a  obtained  from  the  full 

N 

range  of  a  tested.  For  the  test  cases  where  Cjj  was  linear  with  a  (M  ~  3.0), 

these  two  values  should  be  the  same.  These  theoretical  methods  are  usually 
recommended  in  various  handbooks  and  textbooks  because  of  their  success  in 
predicting  C^  at  low  values  of  a.  Their  success  has  been  demonstrated  by 

aW 

several  investigators  at  the  low  values  of  a  and  at  low  supersonic  Mach 
numbers.  The  inadequacy  of  these  theoretical  methods  for  predicting  AC„ 
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without  adding  a  non-linear  term,  such  as  cross-flow  lift,  was  demonstrated  by 
Flax  and  Lawrence3  in  1951.  Cross-flow  lift  for  wings  is  a  concept,  taken 
from  cross-flow  lift  on  cylinders,  which  attempts  to  account  for  the  vortex 
lift.  The  cross-flow  drag  value  used  in  determining  cross-flow  lift  is  basi¬ 
cally  an  experimental  value  obtained  for  a  limited  class  of  wings.3’6  More 
recent  approaches  use  the  concept  of  leading-edge  and  side-edge  suction17’18 
to  account  for  non-linear  lift.  As  far  as  can  be  established  from  the  litera¬ 
ture  this  approach  is  not  applicable  to  the  wing  geometries  of  interest  in 
this  study. 

Returning  to  the  discussion  of  Figure  8,  it  is  noted  that  the  test  values 
of  P  for  thick  housings  is  generally  lower  than  those  for  the  "thin" 

“w 

wings.  A  separate  R.M.S.  curve  for  the  thin  wings  demonstrates  this.  The 
value  of  P  CN  =  4/3  marked  on  the  ordinate  of  Figure  8  will  be  shown  later 

aW 

to  provide  a  reasonable  agreement  with  the  majority  of  test  values  of  AC  for 

W 

the  thick  wings  used  in  this  study,  12  Mach  number -configurational  combinations. 
The  solid  points  shown  in  Figure  8  are  for  test  cases  where  M~  3.0.  In  this 
region  AC^  is  very  non-linear  with  a  at  low  values  of  a.  For  these  cases, 

it  will  be  shown  later  that  Hart's  empirical  method1  provides  good  predictions 
at  the  lower  values  of  a  and  for  some  cases  at  all  values  of  a  tested. 


2.  Thin  Delta  Wings 

A  similar  correlation  plot  of  0  CN  for  the  test  data  for  delta  wings 

“w 

is  given  in  Figure  9  and  is  compared  with  linear  theory  for  these  wings.  In 
this  case  p  C^  is  given  as  a  function  of  P  cot  A.  and  plotted  vs.  1/P  cot. A.  . 

°W 

The  disagreement  with  linear  theory  is  obvious.  Note  specifically 

that  the  test  values  of  P  do  not  tend  to  4  at  P  cot  /!_  =  1  as  predicted 

aW 

by  linear  theory  but  rather  they  tend  to  4  at  P  cot  -A-  =  ®  which  is  in 
agreement  with  predictions  for  rectangular  wings. 

3.  Combined  Correlation  Curve  for  Very-Low-Aspect-Ratio  Wings 

A  comparison  of  the  R.M.S.  curve  of  p  C^  =  F  (1/P  A.R. )  for  thin 

“w 

nearly -rectangular  wings  (Figure  8)  with  the  R.M.S.  curve  of  P  C^  =  F  (1/P  cot-  A 

aW 

for  thin  delta  wings  (Figure  9)  shows  that  the  two  curves  are 

essentially  the  same.  Thus,  one  single  curve  is  proposed  for  predicting 

P  CN  ,  for  both  thin  rectangular  (or  nearly  rectangular)  and  for  delta  wings. 

°W 

The  curve  has  the  same  functional  form  for  P  C^  as  shown  in  Figure  10.  For 

°W 
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the  thick  wings,  P  =  4/3  is 

aW 

not  found  for  thick  surfaces  for 
trend  of  p  for  this  region. 


proposed  for  (1/P  A.R.)  ~  1.5.  Data  were 

the  region  (1/p  A.R.)  ^1.5  to  determine  the 
The  effect  of  wing  thickness  for  ratios, 


t/d,  between  0.2  and  0.1  also  is  not  known;  the  thick  wings  used  in  the 
analyses  had  t/d  ^  0.2;  the  average  "thickness"  for  the  thin  wings  used  was 
t/d  £  0.1. 


In  summary,  the  correlation  curves  of  Figure  10  are  proposed  as  a  simple 
empirical  method  for  obtaining  p  C^  for  very-low-aspect  ratio  wings.  Since 

“w 

in  practice  these  surfaces  are  usually  mounted  on  a  body  of  revolution  the 
mutual  body-wing  interference  should  also  be  accounted  for.  Morikawa's 
factors  are  recommended  for  accounting  for  this  interference  mainly  because 
they  were  used  in  deriving  p  CN  from  test  data.  The  adequacy  of  the  pro- 

QW 

posed  method  for  providing  good  engineering  estimates  of  AC„  =  C  -  CN  at 

W  nbw  nb 

M~  2.5  and  a  to  about  24°  is  demonstrated  in  the  next  section. 


C.  COMPARISON  OF  TEST  VALUES  OF  AC  WITH  EMPIRICAL  PREDICTIONS 
USING  THE  PRESENT  METHOD  ”w 


The  predicted  values 


of  AC 


"w 


are  derived  from  the  empirical  curves  of 


Figure  10  as  follows: 


P  Cx 


w 


57.3  p 


<*W+  V 


(6) 


radian  and  a  is  in  degrees.  These  values  are  compared 

with  the  test  data  obtained  from  CM  -  C„  for  the  29 

nbw  nb 

Mach  number-configurational  combinations  used  in  the  analysis.  Calculated 
values  of  AC  using  Hart's  method1  are  also  shown,  for  the  cases  where  this 
WW 

method  is  applicable,  to  demonstrate  the  adequacy  of  this  method. 

1.  Thick  Wings 

Calculated  values  of  AC„  using  the  present  method,  given  by  Figure  10 

% 

and  Equation  6,  are  compared  with  test  data  from  the  thick  housings  in 
Figures  11  through  15.  Values  of  AC^  obtained  from  Hart's  empirical  corre¬ 
lation  curves1  are  also  shown.  The  comparisons  show,  in  general,  that  the 


where  P  C^  is  per 
.  aW 

in  Figures  11  to  22 
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present  method  with  P  =  4/3  gives  a  good  to  excellent  representation  of 

“w 

the  test  values  of  AC^  to  a  »  24°,  M~  2.5,  for  the  five  housing  configura¬ 
tions  of  Figure  4.  The  predictions  of  the  present  method  are  especially  good 
at  3.0  where  ACN  is  nearly  linear  with  a. 


At  M >C  3.0,  the  data  are  very  non-linear  with  o  at  low  angles-of-attack 
and  Hart's  method  gives  a  better  prediction  than  the  present  method,  see 
Figures  11  and  14.  At  the  higher  Mach  numbers  the  present  method  provides  a 
better  prediction. 

2.  Thin  Nearly-Rectangular  Wings 

The  results  of  the  evaluation  of  the  present  method  for  thin  nearly- 
rectangular  wings  are  given  in  Figures  16  through  19.  The  test  data  are  for 
the  configurations  of  Figure  5.  The  comparisons  again  show  that  the  present 
method  provides  good  predictions.  Hart's  method  also  gives  good  predictions 
in  the  region  of  applicability  of  his  method,  0  A.R.  &  0.8,  but  this  method 
is  not  better  than  the  present  method. 

3.  Thin,  Delta  Wings 

The  present  method  provides  excellent  predictions  for  the  test  data  for 
the  delta  wings  of  Figure  6,  Figures  20,  21  and  22.  Hart's  method  was  not 
derived  for  delta  wings  and  thus  a  comparison  with  this  method  is  not  made 
for  these  wings. 


CONCLUSIONS 


An  empirical  method  is  derived  herein  for  estimating  the  normal  force 
coefficient  (plus  wing-body  carryover),  AC,,  ,  of  nearly-rectangular  thick  and 


thin  wings,  and  of  thin  delta  wings,  of  very-low-aspect  ratio.  The  method, 
in  combination  with  Morikawa's  interference  factors,  gives  good  predictions 
in  the  range  of  Mach  numbers  from  2.5  to  7.7  and  angles  of  attack  to  24°.  For 
near -rectangular  wings  at  M^  3.0,  Hart's  empirical  correlation  curves  of 
Reference  1  are  recommended  for  estimating  ACN  . 
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Fig.  1  Photo  showing  one  type  of  very-low-aspect  ratio,  thick  housings  (with  flow  through) 
used  in  analysis. 
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Fig.  3  Comparison  of  CN  data  from  IRR  E1  housing  with  two  simple  predictive  methods 4>-  0° 
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Fig.  4  Sketches  of  low  A.R.  wings  (housings)  used  in  analysis. 


THE  JOHNS  NOMINE  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL.  MARYLAND 


Ref.  13 


X 

Wx 

(b/2)/d 

A.R.  Sw/S 
(two  panels) 

<b/2)/d 

3 

0.17 

0.08 

1.83 

~T~ 

4 

0.33 

0.15 

3.69 

I 

5 

0.67 

0.31 

7.35 

X/d  =  5.67 


Base 


Wg  (ref.  14):  A.R.  =  0.120;  Sw/S  =  6.285 


0  5.25  12.05  13.8 


Dimensions  in  body  diameters 


Fig.  5  Thin,  low-aspect-ratio,  nearly-rectanyular  wings. 
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Fig.  6  Thin,  low-aspect-ratio,  delta  wings. 
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Fig.  7  Procedure  used  for  linearization  of  test  data  to  derive 
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Fig.  18  Comparison  of  test  and  predicted  values  of  ACN^  for  Wg,  M  -  2.5,  0  =  0,  0A.R.  -  0.275. 
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